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"Abstract 

,We investigate the effect of a uniform background magnetic field on the chiral symmetry breaking in SU (2) Yang-Mills theory 
on the lattice. We observe that the chiral condensate grows linearly with the field strength B up to yfeB - 3 GeV as predicted by 
chiral perturbation theory for full QCD. As the temperature increases the coefficient in front of the linear term gets smaller. In the 
magnetic field near-zero eigenmodes of the Dirac operator tend to have more regular structure with larger (compared to zero-field 
case) Hausdorff dimensionality. We suggest that the delocalization of near-zero eigenmodes plays a crucial role in the enhancement 
of the chiral symmetry breaking. 
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, Electromagnetic interactions are most commonly used to 
probe various properties of other fundamental interactions, in 
particular, of the strong interactions which bind quarks into 
hadrons. In many physically relevant situations, which can be 
pncountered, for example, in astrophysics or nuclear colhsions, 
a reasonable approximation is that of a constant background 
electric or magnetic field. For instance, very strong magnetic 
fields (B ~ 10"' T, ^/eB ~ 1 GeV) could have been produced 
in the early Universe after the electroweak phase transition ^ . 
At the present age of the Universe, the strongest magnetic fields 
can probably be generated in compact dense stars, such as mag- 
netars (B ~ 10'°T, yfeB ~ 1 MeV) [2]. A controllable way to 
subject hadronic matter to a strong electromagnetic field may 
be provided by heavy-ion collisions in very strong laser pulses. 
The latter possibility will be soon realized in such experiments 
as PHELIX at GSI |3] and ELI |4], with estimated values of 
yfeB of order 0.01 MeV (B ~ 10^ T at the intensity of laser ra- 
diation / ~ 10^^ W/cm^). Thus experimentally available field 
strengths are still much smaller than hadronic scale. 

However, very strong electromagnetic fields with strength 
comparable to hadronic scale (B ~ 10'^ T, VeB ~ 300 MeV) 
can be produced in noncentral heavy-ion collisions [5, 6, il. In 
this case, the magnetic field perpendicular to the reaction plane 
is due to the relative motion of the ions themselves. Such a 
strong magnetic field breaks P and CP symmetries and leads to 
charge separation, which manifests itself in a specific correla- 
tions between emitted positively and negatively charged parti- 
cles [6, 7]. Magnetic fields of the same order should also mod- 
ify significantly the properties of QCD vacuum and can even 
change the order of the phase transition to quark-gluon plasma 



i^. In particular, magnetic field couples differently to parti- 
cles with different chiralities enhancing the spontaneous break- 
ing of the chiral symmetry of QCD, (i.e., the symmetry with 
respect to the interchange of massless quarks with opposite chi- 
ralities 1 iL Il2l [13tl)- Thisphenomenon is now ac- 
tively discussed in the literature IMIaEH] in connection with 
the heavy-ion experiments at RHIC and LHC. Another inter- 
esting effect of strong external fields is the enhancement of the 
trace anomaly flT] . 

From a theorist's point of view, the effect of a strong back- 
ground magnetic fields on QCD vacuum is a good test for var- 
ious model descriptions of strong interactions. It turns out that 
magnetic field enhances the breaking of the chiral symmetry. A 
commonly used order parameter for the chiral symmetry break- 
ing is the chiral condensate 



S = 



\qq\0), 



(1) 



which vanishes if the chiral symmetry is unbroken. If the 
strength of magnetic field is larger than the pion mass (which 
is zero for massless quarks) but is still much smaller than the 
hadronic scale, one can use the chiral perturbation theory to cal- 
culate the field dependence of the chiral condensate 1 10]. To the 
lowest order in the magnetic field, the chkal condensate S (B) 
rises hnearly with the field strength B: 



(2) 



where F„ x 130 MeV is the pion decay constant. The result 
(|2|i is specific to QCD: for instance, in the Nambu-Jona-Lasinio 
model or in AdS/QCD dual models the chiral condensate rises 
quadratically with the field strength lll2lll3lll5ll . 



Preprint submitted to Physics Letters B 



November 10, 2009 



There are also other effects of the strong magnetic fields on 
the QCD vacuum. In a separate paper Jl^ we study the polar- 
ization of the spins of the virtual chiral quarks by the strong ex- 
ternal field. The effect of polarization is described by the chiral 
magnetic susceptibility, which is related to phenomenologically 
interesting radiative transitions in hadron physics, to lepton pair 
photoproduction through specific chiral-odd coupling of a pho- 
ton to quarks, and to radiative heavy meson decays ITtIi . The 
chiral magnetization of the QCD vacuum includes two effects 
of the strong magnetic field: the paramagnetic polarization of 
the quarks' spins and the enhancement of the chiral condensate. 
Our paper |16] is devoted to the spin effects, while the present 
article reports our study of the chiral enhancement. 

Below we describe our results of numerical studies of chiral 
symmetry breaking in SU (2) lattice Yang-Mills theory exposed 
to a uniform background magnetic field. Since non-Abelian 
gauge fields interact with Abelian ones only via interaction with 
charged quarks, in our quenched simulations we change only 
the Dirac operator leaving the path integral weight for non- 
Abelian gauge fields intact. Note that the chiral condensate in 



into the Dirac operator by substituting su (2)-valued vector po- 
tential Afj with u (2)-valued field: 



quenched theory grows logarithmically with lattice volume |18 1 
unlike the condensate in the theory with dynamical fermions. 
Thus, strictly speaking, in the thermodynamic limit the chiral 
condensate in the quenched theory is infinite. We argue in Sec- 
tion l2.1l that these logarithmic corrections do not affect the slope 
of the linear dependence of the chiral condensate on the mag- 
netic field (IS) in the quenched limit. Thus, our result should 
have a smooth interpolation from the quenched theory towards 
full QCD. 



2. Numerical results 

2.7. Chiral condensate 

In order to calculate the chiral condensate we use the Banks- 



Casher formula lll9ll . which relates the condensate ([T) with the 
density of near-zero eigenvalues of the Dirac operator D = 



Y, = lim lim 



V 



(3) 



where V is the total four-volume of Euclidean space-time. The 
eigenvalues A„, the eigenmodes t//„ and the density of eigenval- 
ues p ( A) of the Dirac operator are defined by 

Dllfn^An^n, p W) = < J] 5 U - ^) ) , 



where in the quenched approximation the averaging is 
performed over the gauge fields with the weight 
exp i^-S YM [^/j]) ^ YM [^//] is Yang-Mills action. In order to 
implement chirally symmetric massless fermions on the lattice, 
we use Neuberger's overlap Dirac operator flB]. Ultraviolet lat- 
tice artifacts are reduced with the help of the tadpole-improved 
Symanzik action for the gluon fields (see, e.g., Eq. (1) in ll2lil ). 
Uniform magnetic field B in the direction /i = 3 is introduced 



Af ^B/2{x26^i-xi6^2). (4) 

The expression (|4]l is valid in the infinite volume. In order 
to combine it with periodic boundary conditions on the lattice 
we have introduced an additional jc-dependent boundary twist 
for fermions ||22|- In ^ finite hypercubic volume with pe- 
riodic boundary conditions the total magnetic flux tr oug h any 
two-dimensional face of the hypecube is quantized [23]. The 
uniform magnetic field takes quantized values either. 



qB = Ink/L^ 
where 

q = l/3e 



ke. 



(5) 



(6) 



is smallest (absolute value of) electric charge of the quark. 

In our zero-temperature simulations we use two lattice vol- 
umes 14^* and 16"^ and two lattice spacings a = 0.103 fm and 
a - 0.089 fm in order to study finite-volume and finite-spacing 
effects. We also make simulations at nonzero temperature, 
T = 0.82rf, taking the asymmetric lattice 16^ x 6. The critical 
temperature in SU(2) gauge theory is = 313.(3) MeV 11241] . 
The parameters of our simulations and minimal nonzero values 
of the magnetic field are summarized in Table [1] 
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Figure 1 : The dependence of the chiral condensate on the strength of external 
magnetic field at dift'erent lattice volumes and lattice spacings. The solid lines 
are the fits by the function 

Our numerical simulations reveal a gradual enhancement of 
the chiral condensate as the external magnetic field increases, 
Fig H] Qualitatively, this numerical result confirms various the- 
oretical estimations flSEIHllillil. 

The enhancement of the condensate can be traced already in 
the behavior of Dirac eigenvalues. In Figs.|2]and[3]we plot the 
dependence of twelve lowest Dirac eigenvalues for two typical 
configurations of non-Abelian gauge fields on the strength of 
the magnetic field. The stronger the field, the smaller eigen- 
values of the Dirac operator and the larger the density of the 
near-zero modes. The later indicates the enhancement of the 
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Table 1: Parameters of simulations and fitting results for the chiral condensate. 



Parameters of simulations 


Fitting 


results 


TIT, 


Ls 


L, 


P a, fm 


Ljfl, fm 


Vad, fm* 


VeB™„, MeV 


Eg'-*, MeV 


Ab, GeV 
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3.2810 0.103 


1.44 


4.3 
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307(15) 


1.36(13) 
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1.53(11) 
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Figure 2: The 12 lowest eigenvalues of the Dirac operator for typical configu- 
rations of non-Abelian gauge fields vs. the strength of external magnetic field 
for a selected configuration with unit topological charge. This configuration 
contains one zero Dirac eigenmode. 
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Figure 3: The same as in Figure|2]but for a configuration with zero topological 
charge (no zero Dirac eigenmodes). 



chiral condensate according to the Banks-Casher relation OJ. 
Note also that the number of zero eigenvalues of the Dirac op- 
erator remain intact in external magnetic field. Fig. |2] Indeed, 
the Atyah-Singer theorem implies that the number of zero Dirac 
eigenmodes is equal to the topological charge of the gauge field 
configuration (we always have either left of right modes). The 
Abelian fields do not change the topology of the gauge fields, 
so that they do not influence the number of the zero modes. 

We find that the dependence of the chiral condensate on the 
magnetic field strength can qualitatively be described by the 
analytical prediction (|2]). To this end we fit our numerical data 
by the linear function 



2(B) = So 



eB 



(7) 



where Eq = 2(0) and are fitting parameters. According to 
Figure [1] the function O reproduces our data for each set of 
lattice parameters quite well. The the best fit results are given 
in the last two columns of Table [1] 

Before proceeding further we would like to mention a few 
important cautionary remarks. Firstly, our weakest magnetic 
fields are still quite strong as one can see from Table [T] Indeed, 
the strength of the minimal magnetic field, (eBmin)'^^, is still 
greater than the scale imposed by the zero-field chiral conden- 
sate, Xq^^. In this regime the prediction (|2]) of Ref. [ 10] should 
not work, in general, as the linear behavior is expected to be 
realized for much weaker fields. Thus, the linear behavior is an 
unexpected result of our numerical simulations. 

Secondly, the zero-temperature dependence of the chiral con- 
densate E(B) on the magnetic fields slightly deviates from a 
linear behavior at our weakest fields, as one can see in Fig- 
ure[T] Despite deviations are within error, the nonlinear features 
at weaker fields are not excluded by our data. 

Thirdly, our largest fields are much stronger than the chiral 
scale, (eBmax)''^ '^q^ ■ ^^^^ foi" strongest magnetic fields 
that we use we do not see powerlike S (B) oc (eB)^^^ enhance- 
ment, which should be valid on dimensional grounds fl^. The 
reason may lie in the fact that we are studying the quenched 
theory, in which the virtual charged pions are absent, and the 
asymptotical behavior eB)^^^ may be invalid until higher ener- 
gies. However, the observed linear enhancement of the chiral 
condensate even in the absence of the pion loops is intriguing 
feature of the non-Abelian gauge theory. 

At zero temperature we used three different sets of lattice 
data to calculate the field dependence (|7]). The data indicates 
that the finite lattice spacing and finite volume corrections to 



3 



our results are quite small in the studied range of parameters. 
The values for the zero field condensate, So, and for the slope 
parameter A^, calculated at different lattice volumes and spac- 
ings are quite close to each other. For further reference we take 
the values at largest lattice with finest lattice spacing at T = 0: 



S^'=[(320 + 5) MeV]^ 



A^'=(1.53 + 0.11)GeV. (8) 



The numerical value ([8]) of the chiral condensate at zero mag- 
netic field and at zero temperature Eg* agree well with other nu- 
merical estimations in quenched SU{2) gauge theory |25]. The 
numerical value of the slope - also given in Eq. ([HI) - is 
relatively close to the T = result of the chiral perturbation 
theory ©, Ref. fH: 



A* = 



AnF„ 



= 1.97 GeV. 



(9) 



The later fact indicates that the quenched vacuum (studied nu- 
merically in this article) approximates the (electro)dynamics 

of QCD with two massless quark flavours (used in theoretical 
I — f 

analysis of Ref. |10]) surprisingly well. Thus, the quenched 
vacuum can qualitatively describe the enhancement of the chi- 
ral symmetry breaking induced by the external magnetic field. 
Quantitatively, the quenching effects amount to be of the order 
to 30%. 

Strictly speaking, the chiral condensate in quenched theory 
grows logarithmically with lattice volume BlSll : 



s(y) = 2o 



V 

1 -H A In — 

Vo 



(10) 



where Vq is a certain characteristic volume. Numerically, the 
logarithmic divergence in Eq. (fTOt is too small to be seen in 
lattice data. Indeed, in our case changing the lattice volume by 
a factor of two changes the chiral condensate only within the 
range of statistical errors (see TablelTJ. However, this fact raises 
a principal question about establishing the correspondence be- 
tween our values of S and A^ and the corresponding param- 
eters in the full theory. Numerically, we see that at our lat- 
tice volumes the chiral condensate is quite close to its value, 
l}^^ = 276(ll)(16)MeVin the full theory 126]. This means that 
the coefficient A in (fTOl i is quite small 1 18] and the volume Vo is 
close to our lattice volumes shown in Table[T] The origin of the 
logarithmic volume divergence in ( fTOl i is the logarithmic diver- 
gence in the one-loop correction to pion self-energy ]10, 18]. 
From the Euler-Heisenberg lagrangian for a pion in external 
electromagnetic field ]10] it is easy to see that the derivative 
of E over B does not contain this divergence. Thus, one can 
expect that the derivative 



dY.{eB) _ So 
d{eB) 



A2 



(11) 



is finite in the infinite-volume hmit. We see that actually our 
quenched result for the slope coefficient A^ is quite close to the 
analytical prediction ^ for the full theory [10.1 . 

The finite temperature reduces both the value of the chiral 
condensate So and increases the value of A^, Table [1] The last 



fact indicates that the linear dependence on the strength of the 
magnetic field gets weaker with increase of temperature. 

Concluding this Section we notice that in Ref. ]32] the effects 
of the external diagonal gluonic fields on the phase diagram of 
QCD were studied. These fields - which are called the Abelian 
chromomagnetic fields - turn out to enhance the chiral conden- 
sate similarly to the external Abelian magnetic field of ordinary 
electromagnetism. Thus, the strong magnetic fields of gluons 
and the strong magnetic fields of photons have a similar influ- 
ence on the condensate. 

2.2. Geometric structure ofDirac eigenmodes 

There are indications that chiral symmetry breaking is 
closely related with the dominance of certain singular field 
configurations in the vacuum of gauge theories, namely, with 
center vortices and Abelian monopoles (for a review of 

monopoles and vortices see, e.g., Ref. ]28]). 

It has been demonstrated that removal of center vortices from 
gauge field configurations destroys both confinement and chiral 
symmetry breaking ]29]. More detailed studies revealed that 
near-zero eigenmodes of Dirac operator, which, according to 
the Banks-Casher formula (O, are responsible for the emer- 
gence of the chiral condensate, are strongly correlated with cen- 
ter vortices lisol] . Similarly to center vortices, near-zero eigen- 
modes of Dirac operator form complex fractal structures with 
Hausdorff dimensionality between 2 and 3 ]30, 31]. Thus, it is 
interesting to study the effect of background magnetic fields on 
the localization of low-lying Dirac eigenmodes. 
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Figure 4: IPR )12t of Dirac eigenmodes witli /!„ < 500 MeV for different 
strengtlis of the background magnetic field. 

A commonly used measure of localization is the Inverse Par- 
ticipation Ratio (IPR): 

I„^vjd\pl(x) with Jd\p„(x)^l. (12) 

Here p„ (x) - tf/l (x) i//„ (x) is the density of the n-th Dirac 
eigenmode and V is the volume of the system. 

The larger is the IPR, the more localized is the eigenmode i/'n. 
For the function localized in a single point, ij/ (x) ~ 5*'*^ (x - xq), 
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the IPR is equal to the total volume of space V, while for a de- 
localized function, il/{x) ~ 1 / Vv, the IPR is equal to unity. 
The dependence of the IPR on the magnetic field strength and 
the eigenvalue A„ is illustrated in Fig. |4] For not very strong 
field strengths, magnetic field decreases IPR and thus delocal- 
izes Dirac eigenmodes. 
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Figure 5: Three-dimensional cuts of the density of a low-lying Dirac eigen- 
mode for the same configuration of non-Abelian gauge fields are plotted at zero 
magnetic field (left) and at VeS = 1.35 GeV (right). 

In order to visualize how the structure of the eigenmodes of 
the Dirac operator is modified in the external magnetic field, in 
Fig.|5]we have shown the three-dimensional cuts of the four- 
dimensional level surfaces of the densities of two near-zero 
Dirac eigenmodes for the same typical configuration of non- 
Abelian gauge fields with background magnetic field equal to 
zero (on the left) and with VeB = 1.35 GeV (on the right). 
While at zero fields the eigenmode indeed has a highly irregu- 
lar structure and stretches through the whole space, at VeB = 
1 .35 GeV it occupies a compact region and has more or less 
smooth level surfaces. However, at zero field the IPR is still 
significantly higher Why, on the one hand, the magnetic field 
seemingly increases the localization, and, on the other hand, de- 
creases the IPR? This apparent paradox can be resolved if one 
remembers that at zero magnetic field near-zero eigenmodes of 
the Dkac operator are localized on lower-dimensional irregular 
manifolds i31il . which are strongly correlated with the world- 
sheets of center vortices |30]. These manifolds percolate and 
thus fill all the physical space, in the sense that any two arbi- 
trarily distant points have a non-zero probability to belong to 
the same vortex worldsheet. The total volume of these mani- 
folds is, however, zero, and their IPR diverges as some power 
of ultraviolet cutoff llsoi 31]. If, on the other hand, the eigen- 



modes are localized in finite four-volumes, their IPR remains 
finite in the continuum limit. 

In a non-Abelian gauge theory only quarks interact with 
background Abelian fields, and in quenched theory the path- 
integral weight for the gauge fields is not modified. In par- 
ticular, background fields do not affect positions of center vor- 
tices. Since the configurations of center vortices are the same 
for both plots in Fig. |5] we have to conclude that in background 
magnetic field near-zero Dirac eigenmodes become less corre- 
lated with center vortices. For sufficiently strong fields these 



modes are no longer associated with lower-dimensional struc- 
tures with typical size set by the ultraviolet cutoff scale, but are 
rather localized on a physical scale determined by both hadronic 
scale and the strength of magnetic field. 

The fact that at some large VeB (» 2 GeV in our case) the 
IPR starts growing with magnetic field can be explained as fol- 
lows: at this magnetic field and at finite lattice spacing the total 
volume of lower-dimensional manifolds on which the eigen- 
modes can be localized exceeds the localization volume set by 
the magnetic field. Then this effect should disappear in the con- 
tinuum limit, fl — > 0. One the other hand the increase of the 
IPR at high magnetic fields may have another simple explana- 
tion. At very strong fields the structure of the quark eigenmodes 
should dominantly be determined by the magnetic field and not 
by the underlying non-Abelian gauge field. As the magnetic 
field increases the quark eigenmodes tend to get localized in 
the transverse (to the field) directions filling out low Landau 
levels. The radius of the Landau levels shrinks as the magnetic 
field gets stronger and this effect may lead to the observed in- 
crease in the IPR. However, this effects operates at very strong 
gauge fields, at which the chiral condensate should presumably 
be determined by the scale of the magnetic field alone [,10i1 . 

The effect of delocalization of near-zero eigenmodes in the 
background magnetic field can also qualitatively explain why 
the chiral condensate increases with the field strength. Weaker 
correlation between center vortices and near-zero Dirac eigen- 
modes implies an increase of the total volume of space where 
these modes can be localized. Literally, high magnetic fields 
tear away the eigenmodes from the center vortices. Thus, near- 
zero modes become more degenerate and the density of near- 
zero eigenvalues p (/I — » 0) increases, implying growth of the 
chiral condensate (|3]l. 

3. Conclusions 

Summarizing, we conclude that our numerical simulations 
show a clear qualitative evidence that background magnetic 
field enhances spontaneous breaking of chiral symmetry in a 
non-Abelian gauge theory. We confirm that the chiral conden- 
sate rises linearly with increase of the field strength (|7]l, in at 
least qualitative agreement with the prediction Q of the chiral 
perturbation theory [10]. The linear dependence of the con- 
densate on the magnetic field is governed by the parameter Kb. 
Surprisingly we have found that at zero temperature the value 
of the parameter A^ of the quenched theory ([8]) is close to the 
analytical prediction of the chiral theory (|9|l- 

We argue that the slope of the linear dependence of the chi- 
ral condensate on the strength of the magnetic field (fTTT l is not 
affected by the logarithmic volume dependence specific to the 
quenched limit PIS']. 

The linear dependence of the chkal condensate on the 
strength of the magnetic field is an unexpected result of our 



In a free theory with background magnetic field near-zero modes can be 
localized in any point of physical space, and the chiral condensate {3) indeed 
diverges in the continuum limit I91 I22I1 . 
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study because our fields are much stronger than the fields used 
in the analytical investigation of Ref. jlOll . Moreover, the en- 
hancement of the chiral condensate in QCD is caused by the 
coupling of the external magnetic field to the charged pion 
loops OlOll . which are absent in the quenched approach used 
in this paper. Our result indicates that non-Abelian gauge fields 
alone play a significant role in the enhancement of the chiral 
condensate at nonzero magnetic fields. The dominance of this 
contribution may be revealed in further lattice simulations with 
dynamical fermions. 

We found that the dependence of the chiral condensate on the 
magnetic field - encoded in the parameter - becomes much 
less pronounced with increase of the temperature. 

We observe that the (de)localization properties of near-zero 
Dirac eigenmodes are crucial for the enhancement of chiral 
symmetry breaking in the external magnetic field. Our results 
suggest that the microscopic mechanism of the enhancement is 
the delocalization of the near-zero eigenmodes off the center 
vortices which increases the degeneracy of near-zero modes. 
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